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ABSTRACT 

We apply Minkowski functionals and various derived measures to decipher 
the morphological properties of large-scale structure seen in simulations of grav- 
itational evolution. Minkowski functionals of isodensity contours serve as tools 
to test global properties of the density field. Furthermore, we identify coherent 
objects at various threshold levels and calculate their partial Minkowski func- 
tionals. We propose a set of two derived dimensionless quantities, planarity and 
filamentarity, which reduce the morphological information in a simple and intu- 
itive way. Several simulations of the gravitational evolution of initial power-law 
spectra provide a framework for systematic tests of our method. 
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1. Introduction 

The existence of tlie large-scale clustering of galaxies had already been well established 
by the early 1970's mainly due to the pioneering work of [Ibtsuji fc Kihara (1969)] and 



Peebles (1974) who showed that the two-point correlation function for galaxies in the Lick 



and Zwicky catalogues was positive and had the power-law form ^(r) oc r^^'^ on scales 
^ 10/i^^Mpc. Their result was later extended to three-dimensional galaxy catalogues as 
well. Although the clustering of galaxies is now a well-known fact, a complete description 
of clustering which includes its geometrical features has so far eluded researchers. This is 
perhaps due to the fact that the galaxy density field which we observe appears to be strongly 
non-Gaussian. A Gaussian random field is uniquely described by its power spectrum P{k), 
or its two-point correlation function ^(r) since ^(r) and P{k) form a Fourier transform pair. 
This is no longer true for a non-Gaussian field for which ^ must be complemented by other 
statistical descriptors which are sensitive to the structure of matter on large scales. In the 
so-called "standard model" of structure formation, an initially Gaussian density distribution 
becomes non-Gaussian due to mode-coupling and the resultant build up of phase correlations 
during the non-linear regime. These phase correlations give rise to the amazing diversity of 
form, which is characteristic of a highly evolved distribution of matter, and is often referred 
to as being cellular, filamentary, sheet-like, network-like, sponge-like, a cosmic web etc. 
Most of these descriptions are based either on a visual appearance of large-scale structure 
or on the presence of features that are absent in the reference Gaussian distribution which, 
by definition, is assumed to be featureless. Gravitational instability, for example, may cause 
CDM-like Gaussian initial perturbations to evolve towards a density field that percolates 
at a higher density threshold, i.e. at a lower filling factor, than a Gaussian field ( [Melott et| 
al. 19831) . Such distributions display greater connectivity and are sometimes referred to as 



being "network-like" ( [Y^ess fc Shandarin 1996| ). 



In order to come to grips with the rich textural possibilities inherent in large-scale 
structure, a number of geometrical indicators of clustering have been proposed in the past 
including minimal spanning trees ( [Barrow et al. 1985 ), the genus curve ( |Gott ill et al. 



1986|) , percolation theory ([Zel'dovich fc Shandarin 1982| ; ^handarin 1983| ) and shape analy- 



sis (pathyaprakash et al. 1998| ; ^ahni fc Coles 1995| and references therein). A major recent 



advance in our understanding of gravitational clustering has been associated with the ap- 
plication of Minkowski functionals (MFs) to cosmology ( [Mecke et al. 199^ . The four MPs 



Vo, . . . , V3 provide an excellent description of the geometrical properties of a collection of 
point objects (galaxies) or, alternatively, of continuous distributions such as density fields in 
large-scale structure or brightness contours in the cosmic microwave background. The scope 
and descriptive power of the MPs is enhanced by the fact that both percolation analysis 
and the genus curve are members of the family. Additionally, as demonstrated by [Sahni et| 



- 3 - 



al. (1998), ratios of MFs provide us with an excellent "shape statistic" with which one may 
attempt to quantify the morphology of large scale structure, including the shapes of individ- 
ual superclusters and voids. Spurred by the success of MFs in quantifying the geometrical 
properties of large-scale structure we apply the MFs to scale-invariant N-body simulations 
of gravitational clustering, with an attempt to probe both the global properties and the 
individual "bits and pieces" that might make up the "cosmic web" ([Bond et al. 1996|) . 



2. Method 
2.1. Minkowski functionals 

Minkowski functionals (named after Minkowski 1903 ) were introduced into cosmology by 



Mecke et al. (1994]| employing the generalized Boolean grain model. This model associates 



a body with a point process (in our case given by the location of galaxies or clusters) by 
decorating each point with a ball of radius r. The union set of the covering balls is then 
studied morphologically, whereby the radius of the ball serves as a diagnostic parameter 
probing the spatial scale of the body. 

In this paper we shall use the excursion set approach which is applicable to continuous 
fields (that may be constructed from point processes). [Schmalzing fc Buchert (1997)| pointed 
out how to apply Minkowski functionals to isocontours of continuous fields, where the con- 
tour level (the threshold) is employed as diagnostic parameter. The excursion set approach 
inherits two diagnostic parameters, because we may also vary the smoothing scale used in 
constructing the continuous field. 

In three dimensions there exist four Minkowski functionals V^, fi = 0,1,2,3. They 
provide a complete and unique description of a pattern's global morphology in the sense 
of Hadwiger's theorem ([Hadwiger 1957|) . While reducing the information contained in the 



full hierarchy of correlation functions, this small set of numbers incorporates correlations of 
arbitrary order and therefore provides a complementary look on large-scale structure. The 
geometric interpretations of all Minkowski functionals in three dimensions are summarized 
in Table |I[ 

We calculate global Minkowski functionals of the isodensity contours of the density field 
as described in Appendix 0. Furthermore, we separately calculate the partial Minkowski 
functionals of each isolated part of the isodensity contour. Since the total isodensity contour 
is the union of all its parts, the global functionals are given as sums of the partial functionals 
at the same threshold. This follows the spirit of [Mecke et al. (1994)| , where partial Minkowski 



functionals are introduced to measure local morphology for the generalized Boolean grain 
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model. 

Partial Minkowski functionals (PMF) offer the possibility of probing the morphology of 
individual objects, or the object's environmental morphological properties, respectively. We 
expect that this concept will be more powerful when applied to continuous fields at high 
spatial resolution so that the details of structures are not smoothed out. Their application 
to point processes also delivers more direct information. PMF provide a wide range of 
possibilities for morphological studies which we shall explore in a forthcoming paper. 



2.2. Shapefinders 

One important task of morphological statistics consists in quantifying strongly non- 
Gaussian features such as filaments and pancakes. Given the four Minkowski functionals, we 
aim at reducing their morphological information content to two measures of planarity and 
filamentarity, respectively, as has been done for example with various geometrical quantities 
( [Mo fc Buchert 1990| ), moments of inertia ( [Babul fc Starkman 199^ , and cumulants of 
counts-in-cells ( [Luo fc Vishniac 1995[ ). 



Recently, [Sahni et al. (1998)[ proposed a set of shapefinders derived from Minkowski 
functionals. One starts from the three independent ratios of Minkowski functionals that have 
dimension of length. Requiring that they yield the radius R if applied to a ball, we define 

Thickness T := Width W := — ^, Length L := — ^. (1) 

2Vi 7rV2 4V3 



By the isoperimetric inequalities ( pSf ), we have L>W>T for any convex body. 

Going one step further, ^ahni et al. (199"8]| also define dimensionless shapefinders by 



W -T L-W 
Planarity V := — — —, Filamentarity JF := — . (2) 

W+T L+W ^ ' 



Some examples are in order. 



2.3. Simple examples 



Let us consider some simple families of convex bodies in three-dimensional space that 
can take both filamentary and planar shape. 
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A spheroid with two axes of length r and one axis of length Ar has Minkowski functionals 



Vo = yr3A, V, = -r^l + f{l/X)) , V, = -r {X + f{X)) , V, = 1, (3) 



where0 



, arccosx , , 



By varying the parameter A from zero to infinity, we can change the morphology of the 
spheroid from a filament to a pancake via a spherical cluster. 

A different way of deforming a filament into a pancake goes via generic triaxial ellipsoids; 
thereby one of the smaller axes of a strongly prolate spheroid is increased in size until it 
matches the larger axis and an oblate spheroid has been reached. An integral expression can 
be found in (|Sahni et al. 1998|) . 

Yet another transition from prolate to oblate shape is provided by cylinders of radius r 
and height Ar. Here, the Minkowski functionals are given by 

Vo = 7rr^X, Vi = V(l + A), ^ = ^r(7r + A), ^3 = 1- (5) 

A Blaschke diagram, that is a plot of planarity V and filamentarity !F, summarizing the 
morphological properties of these simple convex bodies is shown in Figure |l]. 



3. A set of iV— body simulations 
3.1. Description 

We start from a family of initial power law spectra P{k) oc A;" with n G {—2, —1, 0, +1} 
set before an Einstein-de Sitter background (fi=l, A=0). We conduct numerical experiments 
using a PM code (consult [Melott fc Shandarin 199^ for details). Four sets of phases were 
used for each model, making a total of 16 simulation runs. Each run consists of 128^ particles 
sampled at an epoch well in the non-linear regime. This epoch is chosen such that the scale 
of non-linearity fcni, defined in terms of the evolved spectrum 

= d'kPik) = 1, (6) 

^0 



^Note that for arguments x > 1, one can use the relation iarccosa; — ln[x + y/x^ — 1) to recover an 
exphcitly real-valued expression. 
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is equal to eight in units of the fundamental mode of the simulation box. By using the stage 
fcjii = 8, we make sure that structure is already sufficiently developed on scales much larger 
than the simulation's resolution, while it is not yet influenced by boundary effects. 

Using a cloud-in-cell kernel, these particles were put onto a 256^ grid, which is the max- 
imum value an ordinary workstation can tackle with acceptable time and memory consump- 
tion. Subsequently, the density field was smoothed with a Gaussian kernel oc exp (— a;^/2A^), 
where x is the distance in mesh units and the width A is set to 3. Tests have shown that this 
value both leads to a reasonably smooth field, and preserves at least some detail on smaller 
scales. Throughout the article, we re-scale the density to the density contrast 6, ranging 
from —1 to infinity with zero mean. 



3.2. The global field 

The global Minkowski functionals calculated from the density fields described in the 
previous section are shown in Figure 0. The four different line styles correspond to the 
different spectral indices. Figure ^ shows the Minkowski functionals for the same set of 
models, but instead of the density threshold 6, the rescaled threshold u is used as the x- 
axis, u is calculated from the volume Minkowski functional Vq, that is the filling factor /, 
as described by |Gott III et al. (1987)| . Essentially, its use forces exact Gaussian behavior of 
the volume Vq, by the implicit connection 

4 4* (^) ■ 

Thus the deviations from Gaussianity that are due to changes in the one-point probabil- 
ity distribution function are removed, and deviations due to higher-order correlations are 
emphasized. 

Obviously, the global functionals clearly discriminate between the various models. How- 
ever, in order to make this statement more quantitative, let us take a closer look at the 
individual coherent objects composing the isodensity contours. 



3.3. The largest objects 

At intermediate thresholds, the excursion sets consist of numerous isolated objects. We 
identify them by grouping adjacent occupied grid cells into one object, where adjacent means 
that the cells have a common face. Since the Minkowski functionals of the global field are 
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calculated by integrating over quantities that can be assigned to individual grid cells, the 
partial Minkowski functionals of each object can be obtained at no extra cost once the cells 
belonging to each object have been identified. 

Several plots in Figures and illustrate the behavior of the Minkowski functionals 

of these objects. 

Obviously, the contribution of smaller objects to the volume is almost negligible com- 
pared to the largest one. Note that in all figures, the mean and standard deviation over 
all four realizations are shown instead of the individual curves. It is worth noting that the 
variance is largest in the n = —2 and —1 models, which are dominated by structures on large 
scales and hence show the strongest sample variance. 

The models with various initial spectral indices n show qualitatively a similar behav- 
ior. At small filling factors (high density thresholds) the two largest clusters give negligible 
contribution to each of the global characteristics. Then, at percolation transition the largest 
cluster quickly becomes the dominant structure in terms of volume, area, and integrated 
mean curvature. The second largest cluster also grows at the percolation threshold but just 
a little and then quickly diminishes. The percolation transition is clearly marked in all three 
characteristics of the largest cluster by their sudden growth. However, this transition does 
not happen at a well-defined threshold. Instead, clusters gradually merge into the largest 
objects as the threshold is decreased (the filling factor grows). This continuous transition 
has also been observed using percolation analysis, i.e. the zeroth Minkowski functional alone 
( ^handarin fc Zel'dovich 1983| ; [Shandarin fc Yess 199^ ; [Klypin fc Shandarin 1993| ; ^ahni et 



|al. 1997| ) and is explained by the finite size of the sample. 

Nevertheless, the percolation transitions happen within fairly narrow ranges of the fill- 
ing factor that are clearly distinct for different models in question: The filling factors are 
approximately 0.03 ± 0.01, 0.07 ± 0.015, 0.11 ± 0.015, 0.14 ± 0.015 in the n = -2, -1, 0, +1 
models, respectively. It is remarkable that all percolation transitions occur at smaller filling 
factors than in Gaussian fields (about 0.16) indicating that even in the most hierarchical 
model (n = +1) the structures tend to be more connected than in the "structureless" Gaus- 
sian field. [Pauls fc Melott (1995)| showed positive correlation with networks based on the 



same phases all the way to n = +3. This confirms the conclusion of |Yess fc Shandarm 



(1996)| : The universality of the network structures results from the evolution of Gaussian 



initial conditions through gravitational instability. The Euler characteristic of the largest 
cluster also marks the percolation threshold but in a different manner: before percolation 
it is zero and after percolation it becomes negative in every model, however, in the n = 
and n = +1 models it grows to a small positive peak before becoming negative. All global 
functionals have no particular features at the percolation threshold. 
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3.4. Small objects 

As an example, the Blaschke diagram for the model n = —2 is shown in Figure ||. The 
distributions for the other models look qualitatively very similar and the average quantities 



for other models are presented in Figures § and |10|. Figure § shows that most of the small 
objects are either spherical or slightly planar (two largest dots in Figure H). There is also 
a considerable number of elongated clusters with filamentarities from 0.1 to 0.5. In some 
cases filamentarity reaches large values ~ 1. On the contrary planarity is much weaker, it 
hardly reaches the value of 0.2 (which is partly a consequence of the smoothing). There is a 
hint of a small correlation between filamentarity and planarity: the objects with the largest 
filamentarity also tend to have larger planarity. 

Figures || and |10| display the shapefinders as functions of the cluster mass. The curves 
give averages over the realizations of each model. Apparently, the signal for filamentarity 
is much stronger than for planarity, regardless of the model which is in full agreement with 
Figure The planarity and filamentarity distributions qualitatively look very similar except 
the amplitude. Small objects (5 x 10~^ ^ m ^ 5 x 10""^) display stronger planarity and 
filamentarity for models with more power on large scales. However, for greater masses (m ^ 
5 X 10~^) the situation is reversed: the less large-scale power the greater the filamentarity and 
planarity. If the former seems to be natural and was expected, the latter has been unexpected. 
Both the planarity and filamentarity monotonically grow and reach their maxima at largest 
clusters: Vm ~ 0.1 and Vm ~ 0.5 in all models. As expected the largest objects possess the 
largest planarities and filamentarities, but the independence of the maxima from the model 



again was unexpected. Figures |TT| and show the histograms for the shapefinders V and 
J-', respectively. They clearly show the large difference in the total number of structures: 
the more power on small scales the greater the abundance of clusters. These figures are in 



general agreement with the Figures U and 10 



4. Conclusion and Outlook 

Global Minkowski functionals do discriminate (Figures ^ and Only the n = —2 
model shows slight drawbacks as far as robustness is concerned, but that is due to the 
method's sensitivity to large-scale features of the smoothed density field. The total area, 
mean curvature, and Euler characteristic are sensitive to abundances of the structures and are 
easy to interpret: the more power on small scales (greater n) the more abundant structures 
and therefore the greater the amplitude of the curve. 

Even more valuable information is obtained from looking at the Minkowski functionals 
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of the largest coherent object at each threshold (Figures]^, ||, and 0). All four Minkowski 
functionals of the largest cluster clearly consistently detect the percolation transition. Two 
points are worth stressing: 1) in all models the percolation transition happens at smaller 
filling factors than in the "structureless" Gaussian fields and 2) the more power on the large 
scales (i.e. the smaller n) the smaller the filling factor at percolation. Both conclusions con- 
firm the results of [Vess fc Shandarin (1996)| about the universality of the network structures 



in the power law models with n <1. The results of Pauls fc Melott (1995)] present evidence 
that this should be expected all the way up to n = +3; at n = +4 mode coupling effects 
from smaller scales should begin to fully disrupt the network structure. 

Small objects, on the other hand, give different results. Their abundance discriminates 
well, but is already determined by the difference in the Euler characteristic as well as by 
the total area and mean curvature of the whole contour (Figures ^ and |^) that are also 
sensitive to the abundance of structures. The morphology of small objects as measured by 
shapefinders shows little differences between models so far (Figures |10|, |ll] and |12D. Both 
the maximum average planarity {V ~ 0.1) and the maximum average filamentarity (JF ^ 0.5) 
are reached in the most massive non-percolating objects. None of the model showed ribbon- 
like objects characterized by both large planarity and large filamentarity. We may speculate 
that the smaller objects are ones which formed earlier, are more nonlinear, and therefore 
more decoupled from initial conditions. 

However, all models used Gaussian initial conditions, and evolve under the influence of 
gravity. Hence the similar morphology of the clumps may point towards universal behavior. 
Note that things such as string wakes might produce totally different results. 

The grouping and measurement techniques used in this study may be less accurate 
for small objects than for large clusters. It is worth trying to study the morphology of 
small objects by applying more accurate methods of measuring partial Minkowski functionals 
such as a Boolean grain model ( ^chmalzing 1996 , and a follow-up to this article) or the 



interpolation method of Novikov et al. (1998) generalized to three dimensions. 
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A. Calculating Minkowski functionals of isodensity contours 

Using methods from integral geometry and scale-space theory, |Schmalzing fc Buchert 
(1997)1 developed numerical methods to quantify the morphology of isodensity contours of 
a random field. They present two complementary possibilities for estimating Minkowski 
functionals. 

Given a density field sampled at the grid points of a cubic lattice with lattice constant 
a, Crofton's formula (|Crofton 1868| ) requires counting of the number of grid cells per unit 
volume contained in the excursion set. We obtain 



Vo{v) = 713, (Al) 

VIM =1^(^2-3^3), (A2) 

V2{y) =9^7(^1-2^2 + 3^3), (A3) 

%{^) = ^K-m+r22-n3), (A4) 



where the quantity rij is the number density of j-dimensional elementary cells; to be specific, 
is the number of grid volumes, ^2 counts the grid faces, and rii and denote the numbers 
of grid edges and lattice points, respectively. 

Alternatively, it is sufficient to estimate the derivatives of the random field at the grid 
points and perform a spatial average = ^ / d^x- over Koenderink invariant^ 



V 

=(^Q{u-u))^, (A5) 

={\5{u-v){u^,Yl^)^, (A6) 

w =(^^(^-^) "-"'^:::;r"" )^- (as) 



Careful analysis reveals that both families of estimators are biased by the finite lattice 
constant a encountered with any practical realization of a random field. The deviations are 
of order 0? and can be evaluated analytically for a Gaussian random field (see Winitzki~S^ 
Kosowsky 1997| for some results). It turns out that the two estimates using Crofton's formula 
and Koenderink invariants, respectively, deviate from the true value in opposite directions. 



^The notation is chosen to emphasize the coordinate invariance of the integrand. We use indices following a 
comma to indicate differentiation with respect to the corresponding coordinate, e^fe denotes the components 
of the totally antisymmetric third-rank tensor normalized to €123 — 1. Summation over pairwise indices is 
understood. 
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Hence their difference gives at least a rough idea of the errors associated with binning the 
continuous field onto a finite lattice. 

Yet another method of calculating Minkowski functionals of isocontours was suggested 
by [Novikov et al. (1998)| . However, their interpolation method has so far only been imple- 



mented in two dimensions. 



B. Isoperimetric inequalities and shapefinders 



In his original article, [Minkowski (1903"]] introduced the mixed volumes V{Ki, 
of n convex bodies Ki to Kn- In terms of ordinary volumes of Minkowski sums of the bodies 
we have 



1 " 



(Bl) 



k=l 



where the Minkowski sum is defined as 



ii<...<ik 



K ® L := {x + y\x e K,y e L} . 



(B2) 



In the case of a single body, the mixed volume reduces to its ordinary volume. Furthermore, 
these quantities inherit many interesting properties from the volume and hence play a central 
role in the Brunn-Minkowski theory (see [Schneider 1993| for an introduction). 

A most useful inequality proved by [Alexandrov (1937][ and Penchel (1936)[ states that 



(B3) 



The Minkowski functionals of a body K in d dimensions can be related to mixed volumes of 
two bodies viaP[ 



-V{K,... ,K,B,... ,B), 



(B4) 



^The volume of the j-dimensional unit ball is denoted by luj. Some special values are too = 1, wi = 2, 
0^2 = TT, and uj^ = The general formula reads 
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where B denotes the (i-dimensional unit ball. Consult the book by [Burago fc Zalgaller 



1988)1 for a broad discussion of the Alexandrov-Fenchel inequality and related issues. 



In three dimensions, the Alexandrov-Fenchel inequality (p3| ) leads to two independent 
non-trivial inequalities for the Minkowski functionals, namely 

Vi^ >fK)^2, (B5) 
>i;ViV,. (B6) 

(B7) 



These two inequalities motivate the introduction of the Blaschke diagram ( [Hadwiger 



1955|) . A convex body with Minkowski functionals Vj is mapped to a point {x,y) with 



coordinates 



nVoV2 8V,Vs 



A ball has x = y = 1. For convex bodies, all Minkowski functionals are non-negative, and 
satisfy the inequahties (p5| ), so these points are confined to the unit square. Nevertheless, 
the convex bodies do not fill the whole unit square - the limiting isoperimetric inequality 
has yet to be found ([Schneider 1993| and references therein) . 

The dimensionless shapefinders introduced in equation (0) are related to the isoperi- 
metric ratios via 

l + x' l + y ^ ' 

Obviously, a scatter plot of shapefinders is almost equivalent to the ordinary Blaschke di- 
agram reflected at the point (0.5,0.5), so either method of presentation should convey the 
same morphological information. Throughout this article, we refer to plots of shapefinders 
(P, JF) as Blaschke diagrams. 
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geometric quantity 




V, 




V 


volume 





Vo 


= V 


A 


surface 


1 


V, 




H 


mean curvature 


2 
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Table 1: Minkowski functionals expressed in terms of the corresponding geometric quanti- 
ties. 
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Fig. 1. — The shapefinders of some convex bodies. The sohd hne shows the transition from 
a ball to a filament to a pancake and back to a ball via triaxial ellipsoids of various shapes. 
The dashed line corresponds to cylinders that also undergo the transition from a filament 
to a pancake by varying their height. The dots on the spheroid and cylinder curves indicate 
typical values of the A parameter from Equations and (|^), respectively. 
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Fig. 2. — Minkowski functionals for the evolutionary stage kni = 8 of four different models 
with initial power law spectra P{k) oc /c". The area between two lines gives the one-sigma 
deviations of four realizations from their mean. The solid lines correspond ton — —2, dotted 
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Fig. 3. — Minkowski functionals for the evolutionary stage k^i = 8 of four different models 
with initial power law spectra. Line styles are the same as in the previous figure, but instead 
of the density threshold 6, the rescaled threshold z/ introduced by |Gott 111 et al. (1987)| is 
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Fig. 4. — The following set of figures compares the global Minkowski functionals of the 
whole density field (solid line) to the Minkowski functionals of the largest (dashed) and 
second largest (dotted) objects. The areas indicate the mean and standard deviation over 
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Fig. 5. — The same tendency as before is visible in this case where n = — 1. As before, the 
dominance of large networked structures is clearly visible - for example, the Euler character- 
istic stays well below zero even for high thresholds, which indicates a spongy structure with 
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Fig. 6. — The situation is different in the case n = 0. The largest object sphts up rapidly, 
and is already simply connected at fairly low density contrast {6 0.5). Also, at thresholds 
around percolation, the second largest object becomes comparable to the largest one, so the 
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Fig. 7. — The effect described in the previous figure is even more pronounced for n = 
+1, since this model is dominated by small-scale structures. Contrary to the models with 
negative index n, where the Euler characteristic of the largest object approaches the value 
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Fig. 8. — The shapefinder or Blaschke diagram of the model with n = —2. The area of each 
dot corresponds to the number of coherent objects whose shapefinders (P, JF) he around the 
center of this dot. Note that the shapefinders of all objects at all thresholds are used in the 
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Fig. 9. — The average planarity shapefinder P as a function of the cluster mass. Line styles 
are explained in Figure |[ As in Figure ||, all objects at all thresholds are used to improve 
statistics. 
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Fig. 10. — This plot shows the same type of distribution as the previous Figure but this 
time the average filamentarity shapefinder T is plotted as a function of the cluster mass. 
Again, line styles correspond to models as explained in the caption of Figure 0. 
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Fig. 11. — The distribution of the planarity shapefinder V. The histogram only includes 
objects larger than 10~^ times the simulation box, a value of the order of the smoothing 
volume. Note that although theoretically, the value should lie between and 1, measurement 
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Fig. 12. — This figure displays the same type of distribution as above, but shows the 
filamentarity shapefinder T . Line styles are explained in Figure |^. 



